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ABSTRACT 

We investigate the method for constructing the invariant mass using the M^-assisted on- 
shell (MAOS) approximation to the invisible particle momenta in the cascade decays of a 
new particle resonance produced at hadron colliders. We note that the MAOS reconstruc- 
tion can be defined in several different ways, while keeping the efficiency of approximation 
at a similar level, and one of them provides a unique solution for each event. It is shown 
that the invariant mass distribution constructed with the MAOS momenta exhibits a peak 
at the heavy resonance mass, regardless of the chosen MAOS scheme and the detailed mass 
spectrum of the particles in the cascade. We stress that the MAOS invariant mass can be 
used as a clean signal of new particle resonance produced at the Large Hadron Collider, 
as well as a model-independent method to measure the masses of new particles involved in 
the event. 
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1 Introduction 



The naturalness principle predicts that there is new physics beyond the standard model 
(SM) at TeV scale, which is the scale of the Large Hadron Collider (LHC) experiment 
running now at CERN. It has also been claimed that one of the most important ingredients 
for the new physics model is the existence of a viable dark matter candidate, which is usually 
a weakly interacting massive particle (WIMP). In large classes of the new physics scenarios, 
the WIMP is often stabilized by a discrete parity, which consequently yields a generic 
collider signature with the missing particles always in pairs. Well-known examples such as 
supersymmetric (SUSY), little Higgs, and extra-dimensional models include a WIMP, which 
is the lightest new particle stabilized by a Z2 parity p~H3]. The typical collider signature of 
the new physics model with conserved Z2 parity is the production of a new particle pair, 
decaying to some visible SM particles (multiple leptons and/or jets) plus invisible WIMPs 
in the final state. 

To uncover the underlying physics, one needs the information of the particle properties 
such as mass and spin, which can be revealed from the reconstruction of the signal events. 
However, in hadron colliders, it is generically impossible because the center-of-mass frame 
of the partonic collision is unknown and there are two invisible WIMPs in each event. 
Still, it has been proposed that the new particle masses might be determined by various 
kinematic variables such as the end points of the invariant mass distributions, transverse 
variables, and the techniques using the on-shell mass constraints in various new physics 
processes with the pair-produced new particles [I]. 

In this paper, we consider a new physics event 

p p(p) ->X + U^Y 1 Y 2 + U-> V^xiik) V 2 (q) X 2(l) + U(u), (1) 

where X, Yi, and Xi = 1? 2) are new particles with a priori unknown mass, and U 
stands for visible particles not associated with the decay of X such as jets from the initial 
state radiation. We assume that X has even Z2 parity, while Yi and Xt have odd Z2 
parity. Consequently, X can be resonantly produced and decay into a parity-odd particle 
pair. Vi's denote visible particles and Xi s invisible ones, which yield the missing transverse 
momentum of the event, 

pV = -Pt - - u T = k T + It- (2) 

A similar event topology has recently received a lot of attention. It describes the dilepton 
channel of the top-pair produced in the s-channel mediation of the color sextet bosons [B], or 
axigluon [7] at hadron colliders. In this case, the eight unknown components of two neutrino 
four-momenta can be solved from the six on-shell mass conditions on the top quarks, 

*The event topology is the same as the "antler" diagram, which was studied in [5]. However, here we 
assume that mx is unknown a priori. 
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W bosons, and neutrinos in addition to the two constraints from the missing transverse 
momentum measurements, up to four possible solutions, as well as two combinations due 
to the charge ambiguity on b quarks [8]. In fact, the event topology ([I]) is one of the typical 
signatures of the new physics model with the WIMP stabilized by Z 2 parity, for instance, 
the heavy neutral SUSY Higgs boson (H/A) decaying to a pair of neutralinos (%°) [9] [TO] 
or the n = 2 Kaluza-Klein (KK) Z-boson (Z^) decaying to a pair of n — 1 KK leptons 
(7^) [IT] , producing the final state of several visible SM particles + WIMPs. The masses 
of the new particles involved in the decay chain are generically unknown, and the number 
of decays in the chain might be too short to constrain the unknown masses and invisible 
momenta in many new physics scenarios unlike the aforementioned top-pair process. 

It has been claimed in the literature that it might be possible to measure the particle 
masses by constructing the transverse mass variables even if there are several invisible 
particles in the final state [121411] . However, we note that the end-point position of the 
transverse mass distribution depends on the existence of kinematic configurations, which 
might be forbidden in some models. On the other hand, it has been recently found that 
the Mr2-assisted on-shell (MAOS) method [15], which was introduced to approximate the 
invisible particle momenta in the physics processes with conserved Z 2 parity, can be used 
to measure the SM Higgs boson mass in the dileptonic WW process [16]. We show that 
the MAOS method can be adopted to construct the invariant mass in the case of the event 
topology (jTJ) , thus enabling one to measure the resonance mass in a model- independent 
way. 

This paper is organized as follows. In Sec. EJ we discuss the features of the transverse 
mass variable for the full system, focusing on the behavior of its end point. Then, in 
Sec. El we describe the definition and the types of solutions for the Mt2 variable, which 
is an integral part of constructing the MAOS momenta, as well as known to be useful for 
measuring the masses of the intermediate on-shell states and the WIMP. For each type of 
the Mj-2 solution, the MAOS momenta have distinctive features, which lead to proposing 
various schemes of the MAOS reconstruction, as discussed in Sec. HI where the construction 
of the invariant mass using the MAOS momenta is also described. The comparison between 
the MAOS schemes and the features of the MAOS invariant mass are shown by performing 
a Monte Carlo (MC) study in Sec. [5j We summarize our conclusions in Sec. [6j 

2 Transverse mass of the full system: My 

In this section, we discuss some features of the transverse mass for the event topology 
(pjj. The impossibility of constructing the invariant mass of the decay system due to the 
existence of invisible particles in the final state leads to the proposal of a transverse mass 
variable, which does not use the unknown longitudinal component of the momenta [17] . 
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The transverse mass for the event type (OQ) can be written as 



Ml {Y X Y 2 ) = ml iV2 + ml X2 + 2 (7|p^| 2 + m^p^l 2 + ™ 2 XlX2 - vV 2 ■ vf X2 ) , (3) 

where m Q and pj. denote the invariant mass and the transverse momentum, respectively 
of a = V1V2, XiX2- It is then obvious that Mt(Y{Y 2 ) is bounded from above by mx, thus 
making the determination of nix possible if Mt(YiY 2 ) is correctly constructed event by 
event. It is, however, generally impossible to determine m XlX2 event by event, even though 
the lower bound might be deduced from the knowledge of the event topology and m Xi 
values in some cases, e.g., the SM Higgs boson which decays into two leptonically decaying 
W bosons (h — > WW — > lul'u'). If the lower bound of m XlX2 is determined, it must be 
always true that 

M T {Y X Y 2 ) | mxix2=m mm 2 < M T (y 1 y 2 )U xiX2=m t T e 2 < mx , (4) 

since Mt{YiY 2 ) is a monotonically increasing function of m xlX2 . This fact unties us from 
ignorance of the event-by-event values of m xlX2 , and consequently allows us to construct the 
transverse mass, which can be used to extract the information of mx P2HI3]. For the sake 
of discussion from now on, we focus on the symmetric decay chains, i.e. my 1 = my 2 = my 
and m xi = m X2 = m x . However, we notice that the arguments below can be generalized 
to the case of asymmetric decay chains. 

If m x is unknown and there is no viable theoretical assumption on the lower bound 
of m xx , the best choice will be 

M|(0) = M|(Yiy 2 )| mxiX2=0 = m 2 ViV2 + 2\fa\yf\pT + q T | 2 + m 2 ViV2 - 2(p T + q T ) ■ pV, (5) 

then M max (0) = mx if and only if m x l X2 = 0. The maximum of Mr(0) occurs when Xj's 
are moving parallel to each other in the X rest frame, and the invariant mass of the visible 
particles is minimized. Consequently, the analytic expression of the M^ ax (0) is obtained 
as [13] 

M- ax (0) = £| (A + yjA*-4mi (m^) 2 ) , (6) 
where A = m\ — m 2 x + (m™ m ) 2 and m™ m = min{my i; my 2 } for all the events. If m™ m = 0, 

(77i 2 \ 
1 - -T ' ( 7 ) 
m Y ) 

which shows that the Mt(0) distribution can never reach mx if Xi is massive. The Mr(0) 
variable was found to be useful for measuring the SM Higgs boson mass in dileptonic WW 
decay mode [T2] . 
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On the other hand, in many new physics scenarios, m x is likely to be determined by the 
other kinematic variables such as the Mti kink, which will be described in Sec. El or possibly 
the combinations of the invariant mass end points from various other decay processes. In 
such cases, the more plausible choice of the m x l X2 value will be 2m x (= m Xl + m X2 ) as 



i\4(2m x ) = + Am? + 2w|p T + q T | 2 + ml v J\fa\ 2 + 4m 2 - 2 (p T + q T ) • fa. (8 



V1V2 x V I " I "'ViVa 

But still, the knowledge of m x does not guarantee the saturation of the bound on mx- 
If there is no kinematic configuration in which m XlX2 = 2m x is achieved, the Mr(2m x ) 
distribution will not reach mx, because 

M T {Y X Y 2 ) \ mxx =2m x < M T {Y 1 Y 2 ) \ mxx =m^ < M T (Y 1 Y 2 ) \m xx =m*™ < %• (9) 

The condition for the existence of the kinematic configuration saturating the bound of 
Mt(2?ti x ) up to mx depends on the mass pattern of the involved particles in the decay 
channel, and has been derived in 1 1 3 1 : 



my + m 2 — ( 
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m x < * '-. (10) 

m x 

One can see that it may be impossible to satisfy the above condition if X is too heavy, 
compared with the mass scale of its decay products. In case the condition ( flOl is not 
satisfied, mx is above the upper bound of Mr(2m x ), which depends on the possible range 
of m VlV2 : 

M^ ax (2m x ) = 2m x + m^ 2 



2m * + S|l A ^ 



4m 



m x 



^ v /A 2 -4m 2 ,(m^ in ) 2 I . (ir 



In order to check the properties of M-r(2m x ) discussed above, we have generated MC 
event samples of heavy neutral Higgs bosons (H/A) for two SUSY benchmark points. In 
this model, the heavy Higgs boson decays to a pair of next-to-lightest neutralinos (X2)? 
producing the final state of four charged leptons and a lightest neutralino (x?) pair via a 
three-body process, x° — ^ The detailed description of the chosen model is given in 

Sec. El In Fig. [JJ we exhibit the M<r(2m x ) distribution for benchmark points with relatively 
light (left panel) and heavy (right panel) H/A. One finds that the condition fllOp is satisfied 
in the case of a relatively light H/A scenario, such that the M™ ax (2m x ) corresponds to 
tuh/a- However, it is not satisfied in the heavy H/A case, thus M™ ax (2m x ) is lower than 
rniHi a- This observation shows that the method depends highly on the underlying 
model, i.e. the mass splitting of the on-shell states involved in the event. Hence, one needs 
accurate knowledge of my and possibly the range of my as well as m x to determine mj. 
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Figure 1: Comparison plot of the Mr(2m x ) distributions of the SUSY process H/A — > X2X2 ~~ > 
l + l~Xil + l~Xi- The mass parameters are (rn^o, m x o) ~ (110, 61) GeV for both plots, and m H / A = 
252 (433) GeV for the left (right) panel. See Sec. [5] for a detailed description of the model and its 
simulation. 

3 Alternative transverse mass: Mt2 

In a situation where there are two invisible particles in the event, one may exploit the event 
variable Mt2, which was proposed to measure the particle masses in an event topology such 
as Yi Y 2 — > V\{p)xi{k) V2(g)x2(0; "which is, for instance, the typical collider event of pair- 
produced SUSY particles [18]. Not only is the M T2 useful for measuring the new particle 
masses even in the event topology ([TJ), but it is also an integral part of the definition of 
the MAOS momenta, which will be described in the next section. My 2 can also be denned 
in the asymmetric decay event with my 1 7^ n\y 2 and m Xl 7^ m X2 [19]. However, here we 
consider only the event type with the symmetric decay chains as in the previous section. 
The Mt2 variable is denned as 



Mj-2 = min 



max 



{M«(p T , k T , m x ), i\4 2) (q T , 1 T , m x )}] , (12) 



where M T ^ and M T are the transverse masses of the Y\ and Y 2 systems, respectively, 

(a4 1} ) 2 = m\ + m\ + 2 v /|p T | 2 + m 2 Vi y^l 2 + m\ - 2p T ■ k T , 

(M^ 2) ) 2 = m\ + m\ + 2yWl 2 + m$ 2 y/\l T \2 + m\ - 2q T • 1 T . (13) 

Here, rh x , k^, and It are input trial values of the mass and transverse momenta of the 
invisible particles. They are hypothesized values that parametrize our ignorance of the al- 
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k T , 1 T 



Figure 2: A schematic picture of the balanced (left panel) and the unbalanced (right panel) Mt2 
solutions. 

location of missing transverse momenta into the invisible particle momenta in each event, 
in addition to the true value of m x . The momentum constraint in the minimization au- 
tomatically guarantees that the sum of the hypothesized momenta is equal to that of the 
true momenta, which can be determined event by event. When the input trial mass fh x 
equal to the true mass m x , 



provided that the parent particles Yi are on shell. 

The Mt2 solution for the hypothesized momenta can be classified into two configura- 
tions, unbalanced and balanced. The value of Mt2 is given by MP = M T 2) in the balanced 
configuration, while the unbalanced M^ 2 solution is achieved when the condition for the 
balanced configuration, which will be shown shortly, is not valid. The type of the My 2 

(i) 

solution can be deduced from the value of M T at the hypothesized momenta which give 
the global minimum of Mft (i^j). The global minimum of Mjr (M^) corresponds to 
the stationary point for given values of m Vl ijny 2 ) and m x , 



when = ( y m x /mv 1 )pT and It = (m x /mv 2 )qT, respectively. The Mt2 value of the event 
will be obtained from the balanced configuration when 



M™ 2 ax (m x = m x ) = my, 



(14) 




(15) 




(16) 
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is satisfied at the stationary points. [See the left panel of Fig. [2] for illustration.] For events 
in which the condition ( IT6l) is not satisfied, the Mti is given by the larger value between 
global minima of Mj) at the stationary points, 

M™ = max {(M?) , (M?) } = ( m * + * M f|! > M f!' (17) 
T2 u T ^min'V T /mini \ m Va + m x if M? } < Mp . V ; 

On the other hand, the solution for the balanced configuration needs some nontrivial con- 
sideration. The analytic expression of Mt2 was first derived in [20], then further simplified 
in [21], by considering the event types with vanishing upstream transverse momentum 
(UTM), i.e. u T = in It is given by 

2 



rh 2 x + A T + 



\ \ 2A T - m\ - m V2 



where 

A T = ^J\p T \ 2 + m^ iV /|q T | 2 + m 2 V2 + p T • q T , (19) 

which is invariant under the back-to-back transverse boost. This quantity is closely related 
to the Mct in [23], which has been further generalized to the Mqti in the case of two 
invisible particles in the event [21]. 

The investigation of the Mt2 from the above expressions showed that the end-point 
position of the Mt2 distribution, M™ 2 ax , as a function of the rh x has a kink structure at 
fh x = m x if the invariant mass of the visible particles is not fixed, but has a certain range 
in each decay chain [21]. It was also noticed that the kink structure of M™ 2 ax (m x ) will 
appear when there is a sizable amount of UTM in the events [25]. The Mr2-kink method 
makes it possible to measure m x and my simultaneously even if the decay chain is not 
long enough to constrain all the unknowns in the event. This observation suggests that the 
MT 2 -kink method might also be very useful to measure my and m x in the event topology 
([1]), specifically, if Vi is a set of more than two visible particles and/or the full system of 
FiF 2 is boosted by the UTM. 

If rh x = m Vl = m V2 = 0, Eq. f|T8|) becomes even simpler as 

(M T2 ) 2 U x=0 = 2A T = 2 (|pr||qr| + Pt ■ qr) • (20) 

Interestingly, the authors of (32] have recently claimed that the Mt2 method can be applied 
to measure the mass of the SM Higgs boson in dileptonic WW mode, which is also the 
event type ([!]), and used as an event selection cut to efficiently suppress the background 
processes. 



t Recently, the analytic expression of Mt2 has been derived for some special kinematic configurations. 
See Ref. VM. 
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4 Construction of the invariant mass using MAOS 
momenta 



After discussing the transverse mass variables, we deal with the construction method of 
the invariant mass. In Ref. [15], the authors introduced a systematic approximation to the 
invisible momenta by combining the Mt2 solution with the on-shell relations for generic 
events like Y\ Y 2 — > Vi(p)xi(k)V2(q)x2(l)- The MAOS momenta was originally proposed 
to determine the spin of new particles produced at hadron colliders [T5" |l2"6"ll2"T] . but it was 
later realized that it might also be very useful for the mass measurement [T6||28]. 

The definition of MAOS momenta is composed of two parts, the transverse and the 
longitudinal components. The transverse components of the invisible momenta are set by 
the trial momenta which give the value of Mt2, 



kmaos 
T 



imaos 
1 T 



It, 



(21) 



where kr and It are determined once we minimize max{M^\ M^} in (112)) . among all 
possible trial k T and It satisfying k T + \ T = $ T . The longitudinal and energy components 
are then calculated by the on-shell relations for both x% an d Yi, 



(p 
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rh 2 Yv 



(q + l 1 



m 2 Y2 , 



(22) 
(23) 



where m Xi and •my i (i = 1, 2) are input trial masses of the invisible and the parent particles, 
respectively. The longitudinal components of the MAOS momenta are then given by 



^,maos 



eUp) 
i 



Ap L ± Jp 2 L + E 2 {p)JA 2 - EU P )EUk) 



Bq L ± Jql + E*{q)J& - El(q)E 2 (l) 



(24) 



where 



A = \ (m 2 Yi 
E T {p) = 
E T (k) = 



m 2 Xi - m\ 



Pt ■ k 



maos 
T 5 



Prl 2 + m Vi , 

1,maos|2 i ™2 
\ T I _l_m xi' 



B 

E T {q) 
E T (l) 



^X2 



m 



V-2 



q T • 1 



maos 
T 5 



q r 2 + m 



v 2 > 



(25) 



lmaos|2 i ™2 
l T I + m X 2- 



From the above expressions, it is obvious that both fc™ aos and /™ aos are real if and only if 

\A\ > E T (p)E T (k), \B\ > E T (q)E T (l), 

or equivalently 

m Yl > il4 1} , m Y2 > M { t\ (26) 
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where M^p and Mrp are the transverse masses defined in (fl~3l for = k™ aos , It = l maos 
and the trial masses fh Xi . If the decay chain is symmetric, i.e. rh Xl = rh X2 = m x and 
rhy 1 = rhy 2 = rhy, the above condition reduces to 

m y > max { , M T 2) } , (27) 

so that the MAOS momenta are always real if the value of M™ 2 ax (m x ) is chosen as the 
trial mass of the parent particle for a given rh x . Note that the true mass values my and 
m x automatically satisfy the condition (127)1 because M™ 2 ax (m x = m x ) = my. For the 
end-point events of the balanced Mt2, one has MP = il4 2) = my, and thus k ma,os and 
Z maos correspond to the unique solution of the constraints d2j), (122|) . and ( 123]) . In this case, 
since the true invisible momenta also satisfy the same conditions, they must be equal to 
the unique solution when the true mass values m x and my are inserted in Eqs. (I2"2"j) and 
d23]), 

j^maos /^ rue ^maos ^true (28) 

On the other hand, this is not true for the end-point events of the unbalanced Mt2 because 
M T X) ^ il4 2) , which means that only one side of the MAOS momenta corresponds to the 
true invisible momenta. This argument indicates that the accuracy of the MAOS momenta 
can be controlled by imposing a suitable Mt2 cut, which selects the subset of events near 
the Mt2 end point. Even if m x or my were poorly measured, it was shown that the MAOS 
momenta would provide a good approximation to the true invisible momenta by setting 
m x = and my = M™ 2 ax (m x = 0) [15]. One inevitable problem in this definition is that 
there is generically a four-fold ambiguity on the longitudinal and energy components as 
can be seen in Eqs. (1241) . In the absence of any extra constraints or viable assumptions 
on the kinematic structure of the event, it is clear that there should be no preference of 
one solution to the others. We label this scheme of obtaining the solution of the invisible 
momenta as the first kind of MAOS or MAOS1 to distinguish it from the other schemes in 
what follows. 

In addition to the ambiguity of the solutions, a serious problem will arise when one 
or both parent particles Y t are off shell. The on-shell conditions (123 p can be adopted only 
if the parent particles are on shell, such that the mass values are fixed for all the events. 
To make the MAOS method applicable to the situation where the on-shell conditions are 
not valid, one may consider modifying the on-shell relations. One possible scheme is to 
substitute the event variable M^ 2 , instead of the fixed value my, by 

(p + A; maos ) 2 = M| 2 , (q + / maos ) 2 = M| 2 . (29) 

Then, for the events of the balanced M T2 , one has = = M T2 , and thus /c™ aos 

and /™ aos become unique (MAOS2). This is the scheme which was used to measure the SM 
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Higgs boson mass in [IS] , for both vtih > 2mw an d rrih < 2mw casesE On the other hand, 
for the events of the unbalanced Mt2, i.e. MjP ^ \ a two-fold ambiguity still remains 
on one side of the longitudinal components. 

The only possible way to obtain the unique solution of the longitudinal and energy 
components for all the events is to take 

(p + k maas ) 2 = (M?) 2 , (q + / maos ) 2 = (i\4 2) ) 2 . (30) 

Adopting the above new scheme, the longitudinal components are given by 

7,maos _ Er(k) ^ as _ E T {1) , < 

which are uniquely defined in the events for both balanced and unbalanced Mt2 (MAOS3). 
Note that if the true mass m x is chosen as the input, the MAOS momenta will be equal to 
the true invisible momenta for the end-point events of the balanced Mt2 because M^p = 
(m x ) = m Y corresponds to the true parent particle mass. Although it is 
not true for the events of the unbalanced Mr 2 , one side of the MAOS momenta (fc maos if 
M<p = M£f x = m Y > it4 2) ) still gives the true momenta if the end-point events of Mt2 
were selected. Accordingly, one finds that the accuracy of the MAOS momenta can be still 
controlled by the M T2 cut, which might also be useful for suppressing the SM backgrounds 
in the search for a new physics signal at the LHC |29j . 

Applying the MAOS method to the event topology (jTJ), it is possible to construct the 
invariant mass of the YiY 2 system, 

(p + q + & maos + Z maos ) 2 = (m£ aos ) 2 . (32) 

Then, the successful reconstruction of the invisible momenta by the MAOS method ensures 
that the MAOS invariant mass (m5 aos ) distribution has a peak at the true mass mx, which 
will also be confirmed by the numerical simulation in the next section. The distinctive 
feature of the MAOS method is that it is less model dependent than transverse mass 
variables since the latter highly depends on the mass gaps in the model. Even when 
the resonance particle X is too heavy, the peak position of the MAOS invariant mass 
distribution corresponds to mx, which might be beyond the transverse mass distribution. 
The requisites for constructing the MAOS invariant mass are the assumption of the event 
topology (pQ) and the knowledge of m x , which are also essential for constructing Mr(2m x ) in 
Sec. 121 We also note that the MAOS invariant mass is distinguished by the peak structure, 
and it is generally less vulnerable to the background and momentum smearing effects. The 
clear peak structure of the MAOS invariant mass can be viewed as the smoking-gun signal 
of the heavy resonance produced at hadron colliders. 



* Note that there is only balanced configuration for the SM Higgs decay events, h — > WW — > Ivl'v' . 
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5 Monte Carlo study: heavy SUSY Higgs bosons 



In this section, we illustrate the discussion of the previous section by performing a MC 
study for the decay signal of heavy SUSY Higgs bosons H/A into a pair of next-to-lightest 
neutralino, x®, followed by the decay X2 ~~ >" l + l Xi = e i A*)- This process results in a 
four-lepton plus missing transverse energy final state, 

H/A^xlxl^i + n + r + $t (33) 

As a specific example, we examined two benchmark points in the minimal supergrav- 
ity scenario chosen in [10]. The superparticle mass spectrum has been calculated with 
SOFTSUSY [30], and is given by 

• point A: mu 1 a = 252 GeV, m^o = 110 GeV, m^o = 61 GeV; and 

• point B: uih/a = 433 GeV, m^o = 112 GeV, m^o = 62 GeV. 

We used PYTHIA 6.4 to generate the MC events in the LHC beam condition with the 
proton-proton center-of-mass energy of 14 TeV [31]. For the simple illustrative study, we 
have not considered the detector effects such as the momentum smearing effect and the 
identification efficiency of the leptons. The integrated luminosity is assumed to be large 
enough to measure the particle spectra. The main background process from the SM will 
be ZZ* /'j* with four leptons in the final state. This can be eliminated by requiring large 
missing energy and imposing a Z-veto cut, which rejects events of a dilepton pair with 
the invariant mass near mz- The dominant source of the background in the SUSY process 
is the production of leptons from the decays of neutralinos and charginos, produced by 
squarks and gluinos. In this case, however, the leptons are produced in association with 
jets, so that a jet-veto cut should be imposed to suppress this type of background. The 
direct production of a neutralino pair via the Drell-Yan processes could be challenging 
because it has the same final state. See Ref. [TO] and Chap. 11 of [32] for a detailed study 
of the signature including the backgrounds at detector-level. Here, we do not consider the 
background effect and the event selection cuts except the My 2 cut. 

In the above benchmark points, X2 decays to Xi an d two charged leptons via three- 
body process with an off-shell intermediate Z boson or slepton {rni R = 141 and 406 GeV 
for points A and B, respectively), such that < mi+i- < m^o — m^o. The main difference 
between the benchmark points is the mass of H/A, which is relatively light in point A and 
heavier in point B. One can easily find that the condition ( flO]) is not satisfied in point B, 
whereas it is satisfied in point A. This results in the position of M™ ax (2m x ) being lower 
than rriH/A, as can be seen in the right panel of Fig. [TJ 

To see the characteristic feature of the MAOS momenta, in Fig. [3] we show the dif- 
ference between the reconstructed and the true momenta of Xi f° r point A. The left panel 
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Figure 3: The distributions of k^ aos - k^ me and kf aos - k^ ne for the full event set (left panel), 
and top 10% end-point events of Mt2 (right panel) with the input mass rh x = m^o for the model 
point A. For the MAOS1 scheme, the input mass for the parent particle my is set to m^o. 

includes the distributions of the full event set for the process (133]) generated at the LHC 
condition, while the right panel shows the distributions of the top 10% subset events near 
the end point of Mt2- By definition, each MAOS scheme gives the same transverse MAOS 
momenta. For /c™ aos — k^ ue in the MAOS1 and MAOS2 schemes, we construct their distri- 
butions using all the possible solutions in each event. The result shows that the MAOS1 
scheme is slightly better than the others if one considers the full event set. However, all the 
schemes provide a similar performance if one employs a suitable Mt2 cut. This observation 
and the fact that the MAO SI scheme cannot be adopted when one or both parent particles 
are off shell suggest that the MAOS2 or MAOS3 schemes can be used safely without a big 
loss of efficiency and regardless of whether the parent particles are on shell or not. The 
efficiency of the MAOS momenta will also vary as the detail of the matter content and 
the decay process in the model. It depends not only on the mass spectrum, but also on 
the coupling structure of the particles involved in the cascade decay. The mass gap of the 
decaying particles typically affects the size of the momentum of the visible particles, and 
their particular momentum direction might be forbidden by the helicity correlations. These 
consequently give rise to the different shape of the Mt2 distribution, i.e., how populous 
are the events near the end point. In any case, as pointed out in the previous section, the 
accuracy of the MAOS momenta can be controlled by the Mt2 cut. Although the study 
of the detailed model dependence regarding the efficiency of the MAOS reconstruction is 
beyond our scope here, we stress that it is worthwhile to study before the application to 
the experiments. 
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Figure 4: The distributions of the MAOS invariant mass without the Mt2 cut for point A (left 
panel) and point B (right panel) with the input mass rh x = m^o. In each panel, the MAOS1 
(blue dashed), MAOS2 (black solid), and MAOS3 (red dotted) schemes are adopted to obtain the 
solution of the invisible momenta. For the MAOS1 scheme, the input mass for the parent particle 
mass rhy is set to m^o. 

Next, we proceed to construct the invariant mass of H/A using the MAOS momenta 
following Eq. ( 1321) . In Fig. HI we show the MAOS invariant mass distributions for the 
full event set with the input trial mass m x = m x o while varying the MAOS scheme. One 
can see that all the schemes provide a clear peak structure around run/A- This result 
again indicates that the performance of the MAOS momenta does not depend much on the 
adopted scheme. Figure [5] shows the dependency of the MAOS invariant mass on the input 
trial mass of the invisible particle. Although all the distributions have a peak structure for 
an arbitrary choice of the input trial mass, the peak is located at the true resonance mass, 
m H/Ai only when m x = m x o is chosen. In the right panels of Fig. [5j the top 30% near-end- 
point events of Mt2 are used. This shows that the peak structure of the distributions is 
very clear for the subset of the events, while not changing the position of the peak. Our 
simulation proves that this feature does not depend on the MAOS scheme. 

In real situations, one should consider the combinatorial uncertainty regarding the 
assignment of the visible particles to each chain. For points A and B, the assignment is 
uniquely determined if the event includes two different pairs of opposite-sign same-flavor 
leptons. On the other hand, there are two possible combinations when the event includes 
four leptons with the same flavor, i.e. e + e~e + e~/ / u + / u~/i + / u~. One may select only the 
former type of events to reconstruct the invisible momenta while sacrificing the statistics. 
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Figure 5: The distributions of the MAOS invariant mass for point A (upper panels) and point 
B (lower panels). The full event set is used in the left frames, while the top 30% near-end-point 
events of M?2 are used in the right frames. The MAOS3 scheme is adopted to obtain the solution 
of invisible momenta. 



In a well-known model like the SM, one may calculate the likelihood functions of kinematic 
variables based on the model expectations, then select the combination which gives the 
most likely solution [33]. The other straightforward method is to use the value of Mti in 
the event, as done in J6l[3l]. For a given event, one can construct the Mt2 of all possible 
combinations, then select the combination which gives the smallest value of the Mt2- It 
is actually the same as the M^Qen defined in [2D]. This method can be adopted because 
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Figure 6: The distributions of the MAOS invariant mass for point A (left panel) and point B 
(right panel) with top 30% near-end-point events of Mt2- The MAOS3 scheme is adopted to 
construct the invariant masses, and the input trial mass rh x is set to m^o . In each panel, the true 
combination of four leptons is used for the black solid distribution, while the combination with 
the smallest value of Mt2 is used for the red dotted distribution. 

the Mt2 from the correct pair will have the end point that is definitely related to the mass 
spectrum of the involved particles, whereas there is no such structure in the M^ 2 from 
the wrong pair, which results in a broad distribution in general. Our MC study shows 
that the method is useful for both points A and B. It selects the correct combination with 
68% and 97% efficiency for points A and B, respectively. The different level of efficiency is 
caused by the fact that the decay products are relatively energetic in point B because of 
heavier H/A, and thus it makes the My 2 distribution of the wrong pair even broader. In 
Fig. El we present the MAOS invariant mass for the true combination (black solid), and 
the combination of four leptons with the smallest value of Mt2 (red dotted). If there is 
only one combination, i.e. e + e~fi + fi~ in the event, we use the combination in the MAOS 
reconstruction. The result shows that the overall structure of the distributions remains 
unchanged for both points A and B, such that the method of finding the correct pair using 
Mt2 is successful. 

6 Conclusions 

In this paper, we have examined the possibility of measuring the heavy resonance mass 
by constructing the invariant mass using the MAOS momenta. Regarding the MAOS 
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Table 1: The definition and the number of solutions of the MAOS schemes. 



Definition 



Number of solutions 
Events of Mf 2 Events of M^ 2 B 



MA0S1 Eqs. (|22D and fl23]) four-fold 
MA0S2 Eqs. (|22]) and fl2SD unique 
MAOS3 Eqs. fl22) and (JUD unique 



four-fold 
two-fold 
unique 



reconstruction, we found that various schemes can be defined in order to obtain the MAOS 
momenta, in particular, the longitudinal and the energy components. The MAOS schemes 
are classified by the on-shell equations as summarized in Table [U Although the MAOS 
schemes provide different solutions in general, a similar level of efficiency can be obtained in 
the subset of the events near the end point of Mti- We also note that the MAOS1 scheme 
is not applicable when one or both parent particles Yi are off shell, whereas the MAOS2 
and MAOS3 schemes do not have such a limitation. Using the MAOS reconstruction of the 
invisible momenta, one can construct the invariant mass of the full system. We have shown 
that the peak position of the MAOS invariant mass distribution always corresponds to the 
resonance mass, at which the end point of the transverse mass distribution may fail to point. 
This feature may enable one to deduce directly the mass scale of the heavy resonance even 
in the stage of early discovery, and to measure the mass in a model-independent way. We 
also expect that the MAOS invariant mass distribution can be used as a smoking gun signal 
of the heavy resonance through its clear peak structure, which is generally less vulnerable 
to the background and momentum smearing effects. 
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